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Abstract 

We present here a comprehensive model to describe the bottomonium suppression data obtained from the 
CERN Large Hadron Collider (LHC) at center-of-mass energy of ^snn = 2.76 TeV. We employ a quasiparticle 
model (QPM) equation of state for the quark-gluon plasma (QGP) expanding under Bjorken’s scaling law. 
The current model includes the modification of the formation time based on the temperature of the QGP, 
color screening during bottomonium production, gluon induced dissociation and collisional damping due to the 
imaginary part of the potential between the bb pair. We propose a method for determining the temperature- 
dependent formation time of bottomonia using the solution of the time-independent Schrodinger equation and 
compare it with another approach based on time-dependent Schrodinger wave equation simulation. We find 
that these two independent methods based on different axioms give similar results for the formation time. 
Cold nuclear matter effects and feed-down from higher resonance states of T have also been included in the 
present work. The suppression of the bottomonium states at mid rapidity is determined as a function of 
centrality. The results compare closely with the recent centrality-dependent suppression data at the energies 
available at the CERN LHC in the mid rapidity region. 
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I. INTRODUCTION 


Chu and Matsui [lj employed the concept of 
color screening [2] and presented a model to an¬ 
alyze the quarkonia suppression in quark-gluon 
plasma (QGP). Since then various experimen¬ 
tal [^-0] and theoretical [§-18] studies have been 
carried out on charmonium and bottomonium. 
Bottomonium suppression is considered as a more 
preferred probe in the Large Hadron Collider 
(LHC) experiments since its enhancement in the 
QGP medium is expected to be negligible even at 
the energies available at the LHC. The LHC data 
on bottomonium suppression [lj| at mid rapidity 
have been analyzed using a model based on color 
screening and gluonic dissociation along with colli¬ 
sional damping j20|. However, in the above model, 
the formation time of bottomonia was taken inde¬ 
pendent of the temperature. This is not strictly 
valid. At higher temperatures, the QQ bonding 
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is weakened due to Debye color screening, which 
leads to a decrease in the binding energy and hence 
an increase in the formation time. The formation 
time is a key ingredient for the color screening 
model. Its modification can significantly impact 
bottomonia suppression and consequently it is im¬ 
portant to model the modification of formation 
time due to QGP temperature. The incorporation 
of the temperature effect on the formation time to 
model the quarkonia suppression is currently lack¬ 
ing in the literature. Here, we propose a method to 
estimate the formation time at different tempera¬ 
tures employing the solution of a time-independent 
Schrodinger wave equation. We compare the out¬ 
come of this method with those obtained via the 
solution of a time dependent Schrodinger wave 
equation. It is difficult to obtain the initial wave- 
function for the time dependent Schrodinger equa¬ 
tion exactly. We heuristically arrive at the initial 
waveform. The heuristics and the associated un¬ 
certainties are described in detail in Sec. Ill (B). 

It is expected that the gluonic dissociation and 
collisional damping play a prominent role only af¬ 
ter the quarkonia is formed i.e., after the quarko- 
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nia formation time. However, due to increased 
formation time with temperature, the separation 
between b and b is considerable even before T is 
fully formed. With a temperature-dependent for¬ 
mation time, the role of gluon-induced dissocia¬ 
tion and collisional damping in bottomonia sup¬ 
pression before the bottomonium formation time 
may not be neglected. In view of the above facts, 
it becomes necessary to put forward a suppression 
model based on temperature-dependent formation 
time in order to analyze the experimental data. 


Furthermore, the cold nuclear matter (CNM) 
effects on bottomonia can contribute significantly 
at the energies available at the LHC. Shadowing, 
nuclear absorption, and the Cronin effect are the 
main components of CNM effects. There have 
been studies on the initial state collisional and ra¬ 
diative energy loss in CNM, which can affect the 


nuclear modification factor, Raa 21-t23]. How¬ 


ever, we do not include these initial-state energy- 
loss effects in the current work. It has been argued 
that, as energy increases from the Super Proton 
Synchrotron (SPS) to the Large Hadron Collider 
(LHC) via the BNL Relativistic Heavy-ion Col¬ 
lider (RHIC) experiments, the absorption is ex¬ 
pected to become less and less prominent [24j]. The 
quark and antiquark in the bb pair immediately af¬ 
ter being formed are very close to each other and 
hence the pair behaves almost as a color singlet 
and has negligible interaction with the nucleus, 
leading to almost no absorption. With higher 
energies, the bb pair size is expected to be even 
smaller when it crosses the nucleus leading to even 
smaller absorption. Because we describe p t inte¬ 
grated suppression data, the modeling of the p t 
broadening due to the Cronin effect is irrelevant. 
This leaves shadowing as the principal CNM effect 
that needs to be incorporated. 


There are 
able 


various shadowing models avail- 
in the literature. For the purpose of 
modeling the shadowing effect, we explore two dif¬ 
ferent approaches. In the first approach, we em¬ 
ploy the framework developed by Vogt [2J] which 
is based on the shadowing parametrization first de¬ 
rived by Eskola et al. [29j. From now onwards we 
refer to this shadowing parametrization as EPS09. 
In the second approach, which is a new approach 
proposed by us, we use the shadowing obtained by 


EPS09 for the nuclei whose atomic mass is iden¬ 
tical to or at least close to Np “ rt . . For instance, if 
Npart = 32, then each nuclei provides 16 nucleons 
for the reaction. The effective shadowing would 
mainly be a function of these 16 nucleons. Inter¬ 
estingly, we find that the two approaches yield sim¬ 
ilar results. The shadowing effect, being related 
to the parton distribution function in the heavy 
nuclei, is expected to be an initial-state effect and 
should have the same or similar effect independent 
of the final state of heavy quarkonium (IS' or 2 S 
etc). We use the same shadowing effect in all the 
final bottomonium states. 

In the current work, we describe a com¬ 
prehensive model by incorporating the effect of 
temperature on formation time, color screening, 
gluon-induced dissociation, and collisional damp¬ 
ing along with CNM effects. All these effects are 
modeled analytically. The results of the model 
are compared with the recent centrality dependent 
CMS data [19] at mid rapidity obtained from the 
LHC experiments. The organization of the rest 
of the article is as follows. Section II, briefly de¬ 
scribes the color screening, gluonic dissociation, 
and collisional damping that have been used in 
this work. Section III describes how the tempera¬ 
ture dependent formation time has been modeled. 
CNM effects are also described in this section. Sec¬ 
tion IV gives the results and the comparison with 
the CMS data. Finally, we conclude our work in 
Sec. V. 


II. COLOR SCREENING, GLUON DISSO¬ 
CIATION AND COLLISIONAL DAMPING 


A. Color screening 


The color screening model used in the present 
work is based on the pressure profile [30] in the 
transverse plane and the cooling law for pressure 
based on a quasiparticle model (QPM) equation 
of state (EOS) pij for the QGP. It is described 
in detail in Refs. (20l . 3ll ]. The cooling law for 
pressure is given by 


, , . B C D 

P(r,r) = A+ — + - + —, (1) 
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B. Collisional damping 


where A = -ci, B = c 2 c 2 , C = and D = C 3 . 

The constants < 7 , c 2 and c 3 are given by 
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density described in Refs. 20, |31 


The singlet potential used in this work is given 


by 
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The above constants are determined by using dif¬ 
ferent boundary conditions on pressure and energy 
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where 


Writing Eq. (1) at initial time r = Tj and 
screening time r = r s and combining it with the 
pressure profile |3l], we get the following two equa¬ 
tions: 


rrio is the Debye mass given by mo = 
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p(Ti,r) = A + — + - + —2 = p(n,0)h(r), (2) 
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a eff = yf = 0-63, Nf = 3 = number of 
flavors, oE = 0.47, and a = 0.192 GeV 2 . 


p(r s ,r) = A + 
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where Pqgp is the pressure of QGP inside the 
screening region required to dissociate a partic¬ 
ular T state and it is determined by a QPM 
EOS for QGP medium. The above equations are 
solved numerically and we equate the screening 
time t s to the dilated formation time of bottomo- 
nia t F — 77 >(T) to determine the radius of the 


The collisional damping dissociation time con¬ 


stant is given by T damp = J[^ [Im(V)\ V>] dr [20 
with being the bottomonium wave function. 

We solve the Schrodinger equation to get the 
radial wave function for the 15, 2 S and IP states. 


Gluonic dissociation 


31] The tf(T) is now de- We m °del the gluonic dissociation cross section 


screening region r s 
pendent on the QGP temperature T. 7 = 
is the Lorentz factor corresponding to the trans¬ 
verse energy Ep and the bottomonium mass Mr- 
The expression for survival probability due to color 
screening can be obtained as: 


as 
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with = 0.59, and can be expressed using 
singlet and octet wave functions as: 


where a = 0.5, and Rp and 4> m ax (which is a func¬ 
tion of p t and r s ) are defined in Refs. ( 20 I, 31 


jql' _ 
J nl ~ 


drrg* nl {r)h qV {r). 


(7) 


The above expression for survival probability is in¬ 
tegrated over the range of p t allowed by the CMS 
experiment [19,] to obtain the pt integrated survival 
probability. 


The octet wave function h qV is obtained by 
solving the Schro ding er equation with potential, 
a eff/(8r) [20., 32, [33[. The cross section is then 
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averaged over a Bose-Einstein distribution func¬ 
tion of gluons at temperature T as: 

r , _ 9d f°° dp g P g (Jdiss,nl(Eg) ^ 

diss,nl~7^J^ eEg/T _ 1 , 

with gd = 16 for gluons. 

The net dissociation constant is given by 


r total Fdamp T Pdiss- (9) 

At the very initial time, the bb pair would be very 
close to each other and hence would behave almost 
as a color singlet. During this time, the gluon- 
induced dissociation and the dissociation due to 
collisional damping may be assumed to be negli¬ 
gible. Consequently, we integrate the dissociation 
due to the above two processes from a non-zero 
initial time, to- The survival probability due to 
gluonic dissociation along with collisional damp¬ 
ing is then given by 

/»oo 

S g = exp(-T tota i) dt. (10) 

Jt 0 

We take the value of t 0 to be 0.5 fm. 


D. Net survival probability 

From Ref. [20:], the net survival probability is 
obtained by combining S c and S g 

Qf C C 

O U c Og. 

The expressions for survival probability after in¬ 
corporating feed-down corrections are expressed 
as: 


Sis = 0.6489 S' ls + 0.1363 S' 1P 
+0.1733 S 2S + 0.0416 S' 2Pl 

S 1P = 0.8450 S' ]P + 0.1508 S' 2S + 

0.0041 S' 2Pl 

S 2S = 0.8780 S' 2S + 0.1220 S' P , 

( 11 ) 

where S' nl is the survival probability of the | nl) 
quarkonium states before feed-down is considered, 
while S n i is the survival probability of the | nl) 
state after feed-down. 


III. TEMPERATURE DEPENDENT FOR¬ 
MATION TIME AND CNM EFFECTS 


A. Temperature dependent formation 
time 


We model the temperature as a function of 


Epart &S 20 


T(t X _ T ( N part (bin) \ 1/3 / t QG p \ 1/3 
[q9p) c \N part (bin 0 ) ) {t qgp J [ > 


In the above equation, T c = 0.170 GeV, £qgp is 
the total QGP lifetime taken as 5 fm, and t qgp is 
the current time of the QGP. The temperature is 
inversely proportional to the cube root of proper 
time. The cube root can be seen from the QPM 
as a limiting case when the Reynolds number R 
1 [20]. To model the temperature dependence 
of the formation time we use the real part of the 
singlet potential given in Eq. (5), namely 


V(r,m D ) = —( y l-e~ mDr )- 
m D 

( e ~m D r 
1710 4 - - - 

Solving the Schrodinger equation — 27 §27 + 
V (r, rriE))ip + = Et(ti, l)ip, gives the energy 

eigenvalues Et{u, l) (■n = principal quantum num¬ 
ber, l = azimuthal quantum number) at temper¬ 
ature T for the bottom quark antiquark system 
with reduced mass /i. We calculate the binding 
energy as: 



EbindiT) = Er(n,l) - V(r = 00 ,m D (T)). (14) 

The formation time of T(IS'), T(1P) and T(2S') 



given formation times are greater than the inverse 
of the binding energies in vacuum. I 11 general, the 
formation time is taken to be of the form .. — , 

with K > 1. One can then determine the forma¬ 
tion time at a temperature T\ as: 

Ebind(T 0 ) 


r /( T i) = t/(T 0 ) 


EbindiTi 


(15) 
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The above equation indicates that if the formation 
time in vacuum (i.e., at temperature To = 0 MeV) 
is known, then the formation time at any other 
temperature T { can be calculated. We solve a time 
independent Schrodinger equation and determine 
the T wavefunction ^(r) for various temperatures 
T. Equation (12) shows that the temperature of 
the QGP decreases as the cube root of proper time. 
The question then arises as to what temperature 
needs to be used to determine the formation time 
at a particular centrality. We suggest to use the 
temperature T such that Tf[T(t qgp )\ = t qgp , where 
t qgp is the time for which the QGP has existed. 
The motivation behind this is as follows: When 
the bb pair is formed, it will keep on expanding 
until the stable state of T is formed. Let at QGP 
time to, the temperature of the QGP be T(f 0 )- 
If to < Tf[T(to)\ (given by Eq. (15)), then the 
T has not yet reached its stable state and will 
continue to expand. The desired formation time 
would then be greater than f 0 . As a consequence 
of the above, a point in time will be reached when 
to = Tf[T(to)\. At this point the bb will be at equi¬ 
librium and the T will just get formed. Beyond 
this time the temperature will decrease and the 
T wavefunction may shrink due to string tension, 
but will continue to remain stable unless disso¬ 
ciated by the sea of gluons present in the QGP 
or by collisions with other particles in the QGP. 
The condition Tf[T(t qgp )] = t qgp may not be valid 
for QQ created much after the formation of the 
QGP. As a corollary, the ability of this tempera¬ 
ture dependent formation time model to explain 
CMS data gives an indication that the T enhance¬ 
ment may be negligible at the energies available at 
the LHC. We solve the equation Tf[T(t qgp )\ = t qgp 
by finding the intersection of the two curves t qgp 
and Tf[T(t qgp )\. Figure 1 shows the two curves for 
T(1A) corresponding to the most central bin. The 
point of intersection gives the formation time as 
2.35 fin for T(1A) for the most central bin. 

B. Formation time based on a time- 
dependent Schrodinger wave equation 

We now compare the above method 
of determining the formation time with 
another one based on the solution of a 



FIG. 1: Formation time given by the intersection of 
the two curves t qgp and Tf[T(t qgp )\ for T(15) for the 
most central bin. 

time-dependent Schrodinger wave equation: 
^(aLw) = [n^ ; Q 1 e*' f/ G)G+i-*i)j ip T (r,to), where 
H(ti) is the Hamiltonian with the potential term 
given in Eq. (13), along with the kinetic energy 
and ~ 2 ~~^ terms. The potential is evaluated at 
temperature T(fQ of the QGP at time t,. To 
determine the formation time from simulation, it 
is required to use the following: 

• a criteria to define the condition that the T 
has formed. 

• An initial condition for the simulation. 

For the purpose of determining the formation time 
from the simulation results, we use the follow¬ 
ing criteria: If the T wavefunction ^(r, f) is nor¬ 
malized such that, f ijj^p(r,t)'ijjT(r,t)dr = 1, 

then the mean distance ( r(t )) between the quark 
and anti-quark can be written as: (r(f)) = 

f 'i/j^ r (r,t)r'ijj T (r,t)dr. (Fig. 2) We define the for¬ 
mation time as the earliest time at which the T 
becomes stable. We say that the T has become 
stable when the wavefunction ceases to expand, 
i.e., when = 0. 

Because the effect of temperature can only 
widen the wavefunction and thus increase the for¬ 
mation time (as compared to formation time in 
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vacuum), the simulation for a particular state of T 
is started at a time equal to the formation time tfo 
for that state in vacuum. At any point of time, the 
T wavefunction influenced by the QGP medium 
would be expected to have much larger separa¬ 
tion between the quark and antiquark, as com¬ 
pared to the quark-antiquark separation within 
the bottomonium wavefunction in vacuum which 
has evolved for the same amount of time. This is 
realized by the non zero value of the temperature 
T mit , at which the initial wavefunction is deter¬ 
mined. In fact, for this reason, there is a need to 
apply a lower bound on the value of T init to a value 
higher than 0 MeV. In view of the above argu¬ 
ments, we heuristically fix the lowest value of Ti n u 
to a nominal value of 30 MeV. For the initial wave- 
function at the start of the simulation, we propose 
to use the time independent wavefunction ^(r) at 
the temperature T init , subject to the lower bound 
of 30 MeV, for which the criteria = 0 

exists. Below the temperature T init , the wavefunc- 
tion keeps on expanding. 

For the T(IS') state, we find that the d< '^ f)> = 
0 criteria exists when we start with r/v(r, to) at any 
small value of temperature T, including T — 0 
MeV. Based on the above arguments, we choose 
Tina = 30 MeV. However, for T(2S') and Y(1P), 
the value of temperature varies from T mit = 160 
MeV to Ti n it = 30 MeV depending upon the 

centrality. 

With the above framework for the time depen¬ 
dent Schrodinger equation, the simulation results 
for the three states of T(IS'), Y(2S') and Y(1P) 
are shown in Figs. 3, 4 and 5, respectively. 

We find that the method based on the time- 
dependent Schrodinger equation follows the trend 
of the formation time as a function of centrality 
as determined by our proposed method, i.e., by 
solving the time-independent Schrodinger equa¬ 
tion. Differences in the two curves could possibly 
be due to the uncertainty in the initial conditions. 
For central collisions, the time interval for which 
the time-dependent Schrodinger equation is simu¬ 
lated is large (since the formation time is large), 
and errors in the initial conditions probably have a 
lesser impact. But for peripheral collisions where 
the time interval of simulation is much smaller, the 
impact of the initial conditions is much larger. In 



FIG. 2: Variation of mean separation between the bb 
pair with time for Y(IS') for the most central bin. 



FIG. 3: Comparison of the Y(IS') formation time ob¬ 
tained using the time-dependent Schrodinger equation 
(TDS), our method based on the time-independent 
Schrodinger equation (TIS) and the vacuum forma¬ 
tion time. 

the results section, we describe the T suppression 
by using both the sets of formation time and we 
show that the difference is almost negligible. At 
the same time, we also show that the suppression 
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on the formation time and using the T-dependent 
formation time is significantly different. 

C. CNM effects 

As mentioned in the introduction, nuclear ab¬ 
sorption is expected to be quite small and the 
Cronin effect irrelevant, the main CNM effect is 
therefore shadowing. Thus, only the shadowing 
effect has been used to determine the contribu¬ 
tion of suppression arising due to the CNM effect. 
Shadowing is an initial-state effect and hence it 
is expected to be similar for all the bottomonium 
states [24j. We now describe the two approaches 
that we have used to model shadowing in the cur¬ 
rent work. 


FIG. 4: Comparison of the T(25) formation time ob¬ 
tained using the time-dependent Schrodinger equation 
(TDS), our method based on the time-independent 
Schrodinger equation (TIS), and the vacuum forma¬ 
tion time. 



1. The Vogt approach 

Vogt [24] has computed the shadowing effect at 
a center-of-mass energy of ^snn = 5.5 TeV. We 
employ the same formulation to compute the shad¬ 
owing effect for various centrality bins at yTF/viv = 
2.76 TeV. We use the EPS09 f29[ parametrization 
to obtain the shadowing S l (A,x,p) for nucleus 
with atomic mass number A, momentum fraction 
x, and scale fi. The superscripts i and j refer to 
the incoming and target nuclei, respectively. 

The spatial variation of shadowing 
S p (A,x, fi,r, z), is taken to be a function of 
the shadowing S l (A,x,p) and the nucleon density 
PA{r,z ) as: 

S l p (A, x, /i, r, z) = 1 + N p (S l (A,x,p) - 1) (16) 

/ dzp A (r,z ) 

/ dzp A (0, z) ’ 

where N p is determined by the following normal¬ 
ization condition: 


FIG. 5: Comparison of the Y(1P) formation time ob¬ 
tained using the time dependent Schrodinger equation 
(TDS), our method based on the time-independent 
Schrodinger equation (TIS), and the vacuum forma¬ 
tion time. 


obtained by ignoring the effect of temperature T 


■j j d 2 rdzp A (s)S' r (A,x,ii,r,z) = S'(A,x,p), 

(17) 

with atomic mass number A = 208 for Pb and 
s = \/r 2 + z 2 . The nuclear density p(s) has 
been taken to be the Woods Saxon distribution: 
^( s ) = P° i+ZlV R RA)/d] - The value ^ of A), Ra, d 


and c o have been taken from Ref. 40 for Pb. 
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The suppression factor is defined as the ratio: 


2. Our approach 


RAB^Nparti 


dcr A b! dy 
T AB (b)da pp /dy ’ 


(18) 


where b is the impact parameter and T AB is the 
nuclear overlap function given by 


Tab(V)= / d 2 sdz 1 dz 2 PA(s,Zi)p B (\b- a|,Sa). 

(19) 

In this particular case, both A and B stand for 
atomic mass numbers of Pb. 

From Ref. [4l[, the color evaporation model 
gives 


a AB = I dzi dz 2 d 2 r dxi dx 2 [f g (A,xi,p,r,Zi) 

fg(B, x 2 , p,b-r, z 2 )a gg _ QQ (x 1 ,x 2 , p)}, 

( 20 ) 

a pp = I dx i dx 2 [f g (p, xx,p) 

fgijPi *^2) d)®gg-Q£}A 1 1 *^2) d)\ ■ 

The above formalism excludes the explicit model¬ 
ing of spin and color of the initial and final states. 

In the above expressions, Xi and x 2 are the 
momentum fraction of the gluons in the two Pb 
nuclei at yj snn = 2.76 TeV and are related to the 
rapidity y as: X\ = -J=L=e y and x 2 = J^±- e~ y , 

with M t = \JMy + ( pt) 2 , where Mr is the bot- 
tomonium mass and (p t ) is the mean transverse 
momentum of the bottomonium. 

The function f g (A,x, n,r, z), is determined 
from the gluon distribution function for proton 
f g (p,x,p), by using the following relation: 


• fg(Axi,d,r,z) 

p A (s)S i (A , x u d, r, z)f g (p, xi,p) 

• f’(B,x 2 ,n,b- r,z ) = p B (s)S :t (B,x 2 ,p,b- 
r,z)f g (p, x 2 ,n) 

By using the above expressions, the value of drT fPg B 
and can be computed numerically and, finally, 
R ab can be obtained. 

The value of the gluon distribution function 
f g (p,x,p ) in a proton (indicated by label p) has 
been estimated by using CTEQ6 [42[. 


Repeating the example mentioned in the Intro¬ 
duction, if N part = 32, then each nuclei provides 
16 nucleons for the reaction. The effective shad¬ 
owing is taken to be mainly a function of these 16 
nucleons. In general, 


N / 

S'(A = 20S,x,p,b) = S(A = -^i,x,p), (22) 

where b is the impact parameter corresponding to 
the given value of N part and S^hl = Np “ rt , x, p) 
is the value obtained from EPS09 for nuclei 
with atomic mass N part / 2. We use this proce¬ 
dure to determine the shadowing S Z (A, x, p, b) and 
(A, x, p, b) for the various centrality bins. 

Once S l (A,x, p,b) and Si(A,x,p,b ) are deter¬ 
mined, we determine a AB as 


<y AB = I dx\ dx 2 fg(A, xi, p, b) 
* fg (B, *2, Pz b)O mJ}Q (x i 


(23) 


where 

fg(A Xi, p, b ) = S^A, Xi, p, b)fg(p, Xi, p) 

and 


fi(B,x 2 ,p,b) = S j {B,x 2 ,p,b)f g (p,x 2 ,p). 


We then use Eqs. (18) and (19) to deter¬ 
mine R ab . The shadowing data are not available 
for all the atoms, which leads to the possibility 
that, for a particular centrality bin, the corre¬ 
sponding shadowing data for atoms with atomic 
mass ~ N part /2 would not be available. In such 
cases we use linear interpolation of the available 
shadowing data of the atoms with the closest 
value of atomic mass. We justify the use of linear 
interpolation by the fact that the shadowing effect 
is almost a monotonic smooth function of atomic 
mass (see Fig. 6). 

The atomic mass (more appropriately, we call it 
effective atomic mass) used for each value of N part 
is given in Table I. 

The values of the parameters used in our simu¬ 
lations are given in Table II. 
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TABLE II: Values of the input data used in our sim¬ 
ulation 
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Centrality 

bin 

Npart 

Effective atomic mass 

1 

394 

197 

2 

369 

184 

3 

341 

linear interpolation of 108 and 184 

4 

315 

linear interpolation of 108 and 184 

5 

261 

linear interpolation of 108 and 184 

6 

187 

linear interpolation of 64 and 108 

7 

130 

64 

8 

85.5 

40 

9 

53 

27 

10 

30.3 

16 

11 

16 

9 


T property 

T(1S) 

T(2S) 

T(1 P) 

T(2P) 

Mass (GeV) 

9.46 

10.02 

9.99 

10.26 

t f (fm) 

0.76 

1.9 

2.6 

- 

T diss (MeV) 

668 

217 

206 

185* 


1.2 r 


T 

o CMS data IS 


o 1S: T based formation time 

0 

0 IS: const formation time 


□_B- B -B-B-□ 


FIG. 6: Monotonic and smooth nature of shadowing 
versus atomic mass. The shadowing values shown here 
are for momentum fraction = 0.000312. 

TABLE I: Values of the effective atomic mass and cor¬ 
responding N part . 



IV. RESULTS AND DISCUSSIONS 

The survival probability of the Y(IS') and 
T(2S') states due to color screening versus the 
number of participants are depicted in Figs. 7 and 
8 , respectively. Each centrality region has different 
temperature and the formation time of both the 
T states has been modified based on the temper¬ 
ature corresponding to those centrality bins. The 


FIG. 7: CMS data 19J compared with simulation re¬ 
sults with only color screening for T(IS'). 


experimental CMS data in terms of the nuclear 
modification factor Raa, measured as a function 
of the collision centrality [l9j] are also shown in 
both the plots for comparison. The solid lines indi¬ 
cate the suppression for Y(IS') and T(2S) in Figs. 
7 and 8, respectively, due to color screening with 
temperature-dependent formation time. For com¬ 
parison, the dashed lines in both figures indicate 
the suppression for T(IS') and T(25) bottomo- 
nium states with constant formation time. In the 
peripheral region, where the modification in for¬ 
mation time is small, the suppression due to color 
screening remains similar to the suppression due 
to color screening with constant formation time. 
However, for central bins, where the formation 
time modification is significant, the modification 
in suppression due to color screening is also sig¬ 
nificant. Qualitatively, the suppression now de¬ 
creases while traversing from the peripheral region 
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FIG. 8: CMS data [19j] compared with simulation re¬ 
sults with only color screening for Y(2 S). 


to the central region. At a higher temperature, 
the Y-wave-function stability criterion d< 'jl ]> = 0 
is satisfied at a much later point when the sep¬ 
aration between b and b is much larger. Thus, 
despite weaker bonding due to color screening 
and increased distance between the bottom quark- 
antiquark pair, the quarkonia need not dissociate 
as long as the b, b pair separation is less than the 
critical distance < r >, where = 0. This 

leads to decreased suppression in central collisions, 
where due to higher temperatures, the distance 
< r > at which = 0 happens is larger. It is 

clear from the above figures that the temperature- 
dependent formation time model modifies the sup¬ 
pression upto some extent for T(IS') and is very 
significant for T(2S') suppression. The modifica¬ 
tion of formation time, thus becomes more critical 
for Y(2S). It is also apparent from the above plots 
that the color screening model alone is not able 
to explain the experimental data on bottomonium 
suppression for central collisions. 

In Fig. 9, we introduce gluonic dissociation 
and collisional damping in addition to the color 
screening mechanism. The difference between the 
CMS data and our model prediction has clearly 
decreased in comparison to Figs. 7 and 8 . Our 
suppression results clearly overlap with the mea¬ 


FIG. 9: Survival probability of T(IS’) and Y(2S) ver¬ 
sus Npart by incorporating color screening and gluonic 
dissociation with collisional damping. CNM effects 
are not incorporated. Temperature-dependent forma¬ 
tion time using our proposed method is incorporated 
while calculating color screening. 


sured suppression data within error bars. 

In Fig. 10, the bottomonium suppression with 
all the above three effects, including CNM, is 
shown. This uses the first CNM approach, i.e., 
the Vogt approach. Our simulation results are 
in good agreement with the experimental data. 
The maximum and minimum range of CNM ef¬ 
fects for Y(1S) and Y( 2 S) are indicated by the 
shaded region. The calculations are done using 
EPS09 at leading order and with the temperature 
dependent formation time calculated by solving 
the time-independent Schrodinger wave equation. 
Comparison of the central region for T(IS') be¬ 
tween Figs. 9 and 10 indicates that the mean of 
the experimental data falls between the suppres¬ 
sion without CNM effects (shown in Fig. 9), and 
the mean suppression with CNM effects (shown in 
Fig. 10). The band depicting the possible CNM 
uncertainty overlaps with the error bars. This 
may indicate that the CNM may be on the lesser 
side of the CNM uncertainty band. The CNM 
uncertainty band for Y( 2 S) is pretty small, but 
still a lesser value of the CNM effect continues to 
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FIG. 10: Survival probability of Y(IS') and Y(25) 
versus N par t by incorporating color screening, gluonic 
dissociation and collisional damping along with CNM 
effects. The CNM calculation used here is based on 
the Vogt approach. The bands on the top and bottom 
indicate the maximum and minimum range of CNM 
for T(15) and T(2S), respectively. The temperature- 
dependent formation time using our proposed method 
is incorporated while calculating color screening. 

be in reasonable agreement with the experimental 
data. In Fig. 11, the variation of bottomonimn 
suppression with respect to centrality using the 
temperature-dependent formation time obtained 
by simulation of the time dependent Schrodinger 
equation is shown. Comparison of Fig. 10 with 
Fig. 11 shows that two different approaches to 
determine the temperature-dependent formation 
time give similar results. The differences in the 
peripheral region for T(2S') arises due to the dif¬ 
ferences in formation time in the peripheral region 
between the two methods as shown in Figs. 4 and 
5. 

Figure 12 shows the comparison between the 
two CNM approaches for the variation of the fi¬ 
nal suppression of Y(IS') and T(2S') states with 
N par t■ One can see from the above plot that there 
is little difference in the final suppression employ¬ 
ing the two CNM approaches. Shadowing has 
been shown to vary as A 1//3 [43 : ], Even though 



FIG. 11: Survival probability of Y(IS') and T(2S) 
versus N par t by incorporating color screening, gluonic 
dissociation and collisional damping along with CNM 
effects. The formation time has been determined by 
simulation of the time-dependent Schrodinger equa¬ 
tion. 



FIG. 12: Comparison between the two CNM ap¬ 
proaches for T(IS') and T(25). The subscript 1 
and 2 in CNM\ and CN M 2 refer to the two CNM 
approaches. CNM\ is the ’’Vogt approach”, while 
CNM ‘2 is ”our approach”. 


11 














































we use only the neighboring nucleons, this shad¬ 
owing methodology seems to give results upto a 
reasonable approximation. Figure 13 depicts our 
predicted suppression for Y(1 P) versus N part in¬ 
cluding all the effects employed in Figs. 10 and 11. 
The CNM model used here is based on the Vogt 
approach [24]. Finally, Fig. 14 gives the cornpar- 



FIG. 13: Our prediction for the Y(1P) bottomonium 
state using EPS09 and CTEQ6 after including color 
screening, gluon dissociation, collisional damping and 
the CNM effect. 

ison of the variation of overall Y(1P) suppression 
with respect to centrality determined by employ¬ 
ing two CNM approaches. Again, we see that both 
the CNM methods yield almost the same results. 


V. CONCLUSIONS 

In conclusion, we have presented a compre¬ 
hensive model of bottomonium suppression in 
the QGP medium by combining color screen¬ 
ing, gluonic dissociation, and collisional damp¬ 
ing. Temperature-dependent formation time of 
the bottomonium states is used in the current 
work. We have developed a method for estimat¬ 
ing the modification of the formation time of the Y 
with temperature by solving the time-independent 
Schrodinger wave equation. Temperature depen¬ 
dent formation time has also been determined 



FIG. 14: Comparison between the two CNM ap¬ 
proaches for Y(1P). CNMi is the approach used by 
Vogt. CNM 2 is our approach. 


explicitly by simulating the temporal variation 
of the Y wavefunction using the time-dependent 
Schrodinger equation. We find that the two inde¬ 
pendent methods, based on different approaches, 
give comparable results. The modification of for¬ 
mation time due to temperature modifies the Y 
suppression to a considerable extent in the cen¬ 
tral region and plays a crucial role in accurately 
determining the Y suppression. The quasiparti- 
clc model is employed as an EOS for the QGP 
expanding under Bjorken’s scaling law. The shad¬ 
owing (as a CNM) effect has also been calculated 
by using two approaches, namely, the Vogt [24[ 
approach and our approach. The final suppres¬ 
sion of the bottomonium after taking into account 
the CNM effect is calculated as a function of the 
number of participants and the results are com¬ 
pared with the recent CMS data 19| at the en¬ 
ergies available at the LHC in the mid rapidity 
region. Our simulated results compare reasonably 
well with the CMS data. We also see that both 
the approaches for modeling the shadowing effect 
give similar results. 
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